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ABSTRACT 

We carry out a comprehensive joint analysis of high quality HST/ACS and Chan- 
dra measurements of A1689, from which we derive mass, temperature, X-ray emission 
and abundance profiles. The X-ray emission is smooth and symmetric, and the lensing 
mass is centrally concentrated indicating a relaxed cluster. Assuming hydrostatic equi- 
librium we deduce a 3D mass profile that agrees simultaneously with both the lensing 
and X-ray measurements. However, the projected temperature profile predicted with 
this 3D mass profile exceeds the observed temperature by ~ 30% at all radii, a level 
of discrepancy comparable to the level found for other relaxed clusters. This result 
may support recent suggestions from hydrodynamical simulations that denser, more 
X-ray luminous small-scale structure can bias observed temperature measurements 
downward at about the same (~ 30%) level. We determine the gas entropy at 0.1r v ; r 
(where r v ; r is the virial radius) to be ~ 800 keV cm 2 , as expected for a high tem- 
perature cluster, but its profile at > 0.fr V i r has a power-law form with index ~ 0.8, 
considerably shallower than the ~ 1.1 index advocated by theoretical studies and sim- 
ulations. Moreover, if a constant entropy "floor" exists at all, then it is within a small 
region in the inner core, r < 0.02r v i r , in accord with previous theoretical studies of 
massive clusters. 

Key words: clusters: A1689 - clusters: lensing, X-ray - clusters: DM, gas, temper- 
ature, abundance, entropy 



1 INTRODUCTION 

As the largest gravitationally bound systems displaying a 
range of distinct observational phenomena, clusters of galax- 
ies provide information of central importance in cosmology. 
Total gas and dark matter masses and their profiles provide 
significant insight into the formation and evolution of clus- 
ters and the relationship between baryons and dark mat- 
ter. It has become increasingly clear that much more can 
be learned from careful comparisons of cluster observables 
including galaxy motions, gas properties from X-ray and 
Sunyaev-Zel'dovich (SZ) measurements, and mass profiles 
from lensing distortions and multiple images, providing new 
insight and tighter constraints on the dynamical state of a 
cluster and the nature of dark matter. 

Chandra and XMM observatories have provided very 
detailed information on the physical state of the ubiquitous 
X-ray emitting plasma found in clusters. It has become clear 
that, broadly speaking, there are two classes of clusters, 
those showing some evidence of interaction, evidenced by 
complex structures in the gas, and those for which the gas 
emission is symmetric with a smooth radial variation in tem- 
perature and abundance, indicating the gas is probably re- 



laxed and in a hydrostatic equilibrium with the overall grav- 
itational potential. For lensing work the HST/ACS allows 
the inner caustic structure and central mass distribution to 
be examined in detail (Gavazzi et al. 2002; Broadhurst et al. 
2005a; Sharon et al. 2005). The wide field imagers such as 
on the Subaru and CFHT telescopes permit a statistically 
significant detection of weak lensing distortion and magni- 
fication effects on the background galaxies to be traced out 
to the outskirts of the cluster (Gavazzi et al. 2004; Kneib et 
al. 2005; Broadhurst et al. 2005b). 

In the case of interacting clusters it has proved very in- 
teresting to compare their lensing-based mass distribution 
with their disturbed gas distribution. Clear evidence has 
emerged in the most favorable case of 1E0657-56 (the "bul- 
let cluster" ) that two massive clusters have recently collided 
in the plane of the sky with a high relative velocity, leaving 
the gas lying in-between a bimodal distribution of galax- 
ies and dark matter, traced by weak lensing (Markevitch 
et al. 2002; Clowe, Gonzalez & Markevitch 2004; Bradac et 
al. 2006) . This particular case demonstrates that the bulk of 
the mass is dark and relatively collisionless, as anticipated in 
CDM dominated cosmogonies (Markevich et al. 2004, Clowe 
et al. 2006, Milosavljevic et al. 2007, Randall et al. 2007), al- 
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though the estimated relative velocity between the two mas- 
sive components may be exceptional in the context of ACDM 
simulations (Hayashi & White 2006). Other such examples 
are coming to light, with collisions closer to the line of sight 
(Czoske et al. 2002, Jee et al. 2007, Dupke et al. 2007). A de- 
tailed lensing and X-ray study of a larger sample of interact- 
ing clusters by Okabe & Umetsu (2007) spans the full range 
of dynamical interaction, from premergers where the gas is 
clearly unaffected by the mutual gravitational attraction, to 
cases where both mass components are still readily distin- 
guishable but the gas heavily disrupted and shock heated, 
and finally the postmerger phase where the gas shows only 
local signs of interaction, with a relatively small degree of 
substructure visible in the dark matter as probed by lensing. 

Strong lensing based masses of the central regions of 
massive clusters have often been significantly higher than 
central masses deduced from X-ray analysis, by factors of 
~ 2 - 4 (Miralda-Escude & Babul 1995, Wu & Fang 1997; 
Allen 1998, Wu et al. 1998, Voigt & Fabian 2006). Obvious 
reasons for the different masses - in addition to modeling 
and intrinsic observational uncertainties - include possible 
deviations from hydrostatic equilibrium, sphericity and gas 
isothermality (see Allen 1998). Lensing estimates are natu- 
rally biased upward in the case where the gravitational po- 
tential is significantly elliptical with the major axis preferen- 
tially aligned towards the observer. Evaluation of this effect 
in the context of ACDM simulations shows that the inferred 
concentration of the cluster profile, measured in terms of the 
Navarro, Frenk & White (1996, hereafter NFW) model, can 
be enhanced by up to 20% by this bias in the worst cases 
(Hannawi et al. 2005, Oguri et al. 2005), falling well short 
of the reported discrepancies. 

Motivated by the need to improve the precision of mea- 
surements of gas and DM profiles of clusters, we have begun 
a program of joint X-ray, strong and weak lensing analy- 
ses of several clusters for which we can combine high qual- 
ity resolved observations. The advantages of a simultane- 
ous analysis of X-ray and lensing data are clear, given that 
strong lensing measurements yield the total mass profile in 
the inner cluster core while X-ray and weak lensing mea- 
surements cover a much larger region of the cluster. The 
increasing quality and degree of detail of such data allows a 
more model-independent determination of the relevant pro- 
files along with redundancy so that self-consistency can be 
checked. Under the assumption of spherical symmetry and 
hydrostatic equilibrium the projected temperature and gas 
density profiles as well as the total surface mass density pro- 
file may now be derived directly without resorting to as- 
sumed models or simple parameterizations of the profiles. 

Here we apply a model independent approach to derive 
the density profile of a relaxed cluster from a simultaneous fit 
to both the X-ray and lensing data. We apply our technique 
to A1689 (r V i r ~ 2 h _1 Mpc), a rich and moderately-distant 
(z = 0.183) cluster that has been extensively observed in the 
optical, near IR and X-ray regions, and is the first cluster 
in our sample. The cluster has a cD galaxy whose center 
is within ~ 1.5" of the X-ray centroid; this fact, and the 
low degree of X-ray ellipticity (e ~ 0.08, Xue & Wu 2002, 
hereafter XW02) indicate that the cluster is likely to be well 
relaxed and nearly spherical. Previous estimates of its mass 
were obtained from the analysis of observed arcs and arclets 
produced by strong lensing (Broadhurst et al. 2005a), from 



the distortion of the background galaxy luminosity function 
and number density (Taylor et al. 1998, Dye et al. 2001), and 
from weak lensing observations (Broadhurst et al. 2005a,b, 
Medezinski et al. 2007). Here we derive the X-ray surface 
brightness and temperature profiles from an analysis of the 
full set of Chandra observations, which can be compared 
with a smaller subset of Chandra observations analyzed by 
XW02, and an independent study based on XMM (Anders- 
son & Madejski 2004; hereafter AM04). 

The X-ray and lensing observations and data reduction 
are described in Section 2, followed by a detailed account 
of the spectral and spatial data analysis in Section 3. In 
Section 4 we describe the methodology of deriving the gas 
and mass profiles, and in Section 5 we present the results 
of our deduced gas, total mass, and entropy profiles. Our 
results are discussed and assessed in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 X-ray measurements 

A1689 was observed by Chandra during three non- 
consecutive periods. These observations were made mainly 
with the onboard Advanced CCD Imaging Spectrometer in 
2x2 imaging array (ACIS-I) mode. Table[T]gives a summary 
of the data we have analyzed including the Good Time Inter- 
val (GTI, i.e., exposure time after all known corrections were 
applied). We reduced all data using the following release 
of data reduction software: Chandra data analysis software 
package CIAO 3.3, with the updated complement calibra- 
tion database CALDB 3.2. Q All observations were reduced 
from the level 1 stage in order to achieve a better modeling 
of instrument gain and quantum efficiency. The event grades 
which we used are GRADE=0, 2, 3, 4 and 6. Periods of back- 
ground flaring were removed using the CIAO task "lc_clean" . 
We removed bright sources using the tool "wavdetec" with 
the default parameters: scales="2.0 4.0" and sigthresh = 
10" 6 . Results of setting scales="1.0 2.0 4.0 8.0 16.0" were 
checked for Observation ID 540; no other bright sources 
were identified. We tailored the background files to the cor- 
responding data sets. Identical spatial filters were applied 
to the sources and the background data sets. Spectra and 
responses were extracted using the tool "specextract" . All 
channel count rates were combined into one spectrum (using 
the FTOOL task MATHPHA, and matching response ma- 
trix and ancillary response files, using the tools ADDRMF 
and ADDARF, respectively), weighting individual exposures 
by their respective integration times. We then binned the 
combined counts so that there were at least 25 counts per 

1 We followed the threads for data prepara- 
tion "Analysis Guide: ACIS Data Preparation" 
http://cxc.harvard.edu/ciao/guides/acis_data.html, and for 
extended sources "Analysis Guide: Extended Sources", 
http://cxc.harvard.edu/ciao/guides/esa.html Due to the 
risk that some cluster emission extends over the entire image 
we took the background from the ACIS "Blank-Sky" Back- 
ground files compiled by Markevitch (2001), and followed 
the thread "Using the ACIS "Blank-Sky" Background Files", 

| htt p : / / cxc .harvard . edu/ciao/ threads / aeishackground / 1 and 

Maxim's cookbook, http://cxc.harvard.edu/cal/Acis/Cal_prods/ 
bkgrnd/acisbg/COOKBOOK. 



A 1689 Mass and Gas Profiles 3 



bin. The center of the cluster was found by IRAF to be at 
13 h ll m 29. s 575 - 01°20'27."59, in agreement with the posi- 
tion determined by AM04 (13 h ll m 29. s 4 - 01°20'28"). 

The spectrum of each observation was first checked in 
order to identify possible problematic features, and to as- 
sess the consistency between the model parameters obtained 
from the different observations. We reduced the spectrum 
within a circular radius of 3' (corresponding to a physical 
radius of 387/i _1 kpc), and binned each spectrum to have at 
least 20 counts per bin. Using XSPEC we fitted the 0.3 - 10 
keV data to an optically thin thermal plasma model with 
Galactic photoelectric absorption, WABS(MEKAL), with 
Nh = 2 • 10 20 cm -2 , the mean absorption along the line 
of sight to A1689 (Dickey & Lockman 1990). The resulting 
parameter values were consistent between the fits, but the 
fit quality varied, with Xr = 1-231, 1.366, and 1.478 for ob- 
servation ID 540, 1663, and 5004, respectively. Deviations 
between the data and model were high in two energy bands. 
Indeed, in the 7—10 keV band the count rate is significantly 
higher than the model prediction, as noted also by AM04 
(and as also seen in fig. 2 of XW02). This could partly be 
due to a background of high energy particles, but filtering 
the background of observation ID 540 using a smaller time 
binning did not lower the difference in the count rate. Un- 
corrected instrumental effects can also be invoked, such as 
an imprecise correction for the contaminating lines from the 
external calibration sourc^E 

The second problematic band was 0.3 — 0.5 keV, where 
the data in all three observations are higher than the model 
values (even if absorption is ignored). As mentioned in 
AM04, there is extra absorption caused by molecular con- 
tamination of the ACIS optical blocking filters which causes 
the data to be lower than the model. A correction is imple- 
mented in the analysis software (starting with version CIAO 
3.0), but because of the large uncertainty in the ACIS gain 
at energies below 350 eV, the recommended procedure is to 
ignore events in the 0.3-0.35 keV banc0. We do not know if 
a substantial uncertainty extends also to 0.35 — 0.5 keV, so 
to be safe, we ignored the 0.3 — 0.5 keV data, and used only 
the 0.5 — 7 keV measurements. 

The change in energy interval for the ID 5004 observa- 
tion resulted in only < 1% difference in values of the fitted 
parameters. This "band stability" increases the confidence 
in the reliability of the results. However, the temperature is 
somewhat sensitive to the specific value used for the Galactic 
absorption. Fits with Nh as a free parameter yield a reason- 
able value, 1.32 • 10 20 cm -2 for observation 540, but a value 
which is close to zero for the two other observations. Also, 
kT is higher by ~ 1 keV than in the fit with Nh fixed at its 
observed value (2 ■ 10 20 cm -2 ). This is further discussed in 
section 6 below. The measured flux and surface brightness 
are insensitive to this change, since the overall normaliza- 
tion does not change when the absorption is taken to be a 
free parameter. 

More Chandra observations of A1689 have recently be- 
come public, namely observations ID 6930 and 7289, total- 
ing 80 ks. We have checked the quality of these new data by 



2 For more on this, see http://cxc.harvard.edu/cal/Acis/ 
CaLprods /bkgrnd / current / 

3 http://www.astro.psu. edu/users/chartas/xcontdir/xcont.html 



deriving the temperature profile using CIAO 3.4. The de- 
duced value is systematically offset higher by 1.5 — 2 keV 
than the value obtained from the earlier observations, and 
from values obtained from measurements with other instru- 
ments. Also, the profile from the new observations does not 
show a decrease at large radii but rather an increase. The 
problem is even more severe using CIAO 3.3 and is due to 
improper background files. These results indicate that there 
still is a problem with the matched background files. We 
have therefore not included these observations in our work. 

2.2 Lensing measurements 

Analysis of strong lensing measurements of A1689 was car- 
ried out using deep HST/ACS images with a total of 20 or- 
bits shared between the GRIZ passbands. Over 100 lensed 
background images have been identified, corresponding to 30 
multiply-imaged background galaxies, including many radial 
arcs and small de-magnified images inside the radial critical 
curve, close to the center of mass (Broadhurst et al. 2005a). 
For a given lens model Broadhurst et al. projected the lensed 
images onto a sequence of source planes at various distances. 
They then generated model images by lensing these source 
planes and comparing the detailed internal structures of ob- 
served images falling near the predicted model positions, 
where the unknown source distance is a free parameter for 
each source. As new images were identified they were in- 
corporated into the lens model to refine it, enhancing the 
prospects of finding additional lensed images. This relatively 
rich lensing field allows a good mapping of the mass in the 
inner core. 

At larger radius the statistical effects of weak lens- 
ing have been used to explore the entire mass profile of 
A1689 using wide-field V and I-band images taken with 
Subaru/ Suprime-cam (Broadhurst et al. 2005b). In practice 
this work is difficult, requiring careful analysis of large sets 
of wide field images, with corrections for seeing, tracking 
and instrumental distortion. Broadhurst et al. obtained a 
firmer constraint on the mass profile by examining the effect 
of lensing on the background number counts of red galax- 
ies, as advocated by Broadhurst, Taylor & Peacock (1995). 
Broadhurst et al. found a clear detection of both the weak 
distortions and a deficit in the number counts to the limit 
of the Subaru images, which corresponds to an outer radius 
of r ~ 1.5 h _1 Mpc (Broadhurst et al. 2005b). 

Combining the inner mass profile derived from strong 
lensing with the outer mass profile from weak lensing we 
see that for A1689, the projected mass profile continuously 
flattens towards the center like an NFW profile, but with 
a surprisingly steep outer profile (Broadhurst et al. 2005b, 
Medezinski et al. 2007) compared with the much more dif- 
fuse, low concentration halos predicted for massive CDM 
dominated halos, (e.g., Bullock et al. 2002). 



3 SPECTRAL AND SPATIAL DATA ANALYSIS 

The X-ray flux was measured in two bands - 0.5 — 9 keV, 
which was used to check overall consistency with previous 
observations, and the narrower band 0.5 — 7 keV used in 
this work (see section 2.1). In a 3' aperture the measured 
flux in the latter band was Fo.s-rkeV = (2.17 ± 0.01) ■ 
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Table 1. Chandra Observation Log for A1689 



Obs. ID 


Start time 


Mode 


Duration 


Good time interval (GTI) 








[ks] 


[ks] ACIS 


540 


2000-04-15 04:13:33 


FAINT 


10.45 


10.32 


1663 


2001-01-07 08:18:34 


FAINT 


10.87 


10.62 


5004 


2004-02-28 07:18:29 


VFAINT 


20.12 


18.44 








41.44 


39.38 



10 11 erg cm 2 s 1 , slightly lower than in the wider band, 
F .5-9keV = (2.44 ±0.01) ■ 10" 11 erg cm" 2 s -1 . The corre- 
sponding luminosities were Lo.5_7fc e v = (1.002±0.004) • 10 45 
h" 2 erg s~\ and L .5- 9 fc e v = (1-128 ± 0.004) ■ 10 45 h" 2 
erg s~ . These values are consistent with previous find- 
ings. XW02, who analyzed some of the Chandra data, ob- 
tained -Fo.s-iofcev = (2.7 ± 1) ■ 10" 11 erg cm" 2 s" 1 and 
Lo.s-iokev = (1.03 ± 0.38) ■ 10 45 h -2 erg s" 1 . Our flux is 
also consistent with that measured by XMM (AM04) and 
ROSAT (Ebeling et al. 1996). 

The degree of ellipticity of the cluster X-ray emission is 
not only a good indication of its 3D morphology, but also 
gives some indication of its dynamical state. These impor- 
tant properties are of particular relevance to our work which 
is based on the assumptions that the gas is roughly spheri- 
cally symmetric and that the cluster is in a state of hydro- 
static equilibrium. Using the SExtractor utility we estimated 
the ellipticity of the X-ray emission to be e = 0.083 ± 0.002, 
which is in very good agreement with the value reported by 
XW02. 

We measured the cluster gas temperature from the spec- 
tral measurements (Figure [1} • To check consistency of the 
mean temperature T with previously determined values, we 
analyzed the combined dataset in the 0.3 — 9 keV band 
(which is, as noted above, wider than the more uniform 
0.5 — 7 keV dataset used in the rest of our work). The 
fit (with Xr = 1-14) yielded kT = 9.36 ± 0.18 keV, and 
a metal abundance A — 0.4 ± 0.04 in solar units, consistent 
with kT = 9.02t[;- 4 keV (Mushotzky & Scharf 1997), and 
kT = 8.2 - 10 keV, A = 0.2 - 0.49 (XW02). In the 0.5 - 7 
keV band, the corresponding values are kT = 9.35 ± 0.18 
keV, and A = 0.41 ± 0.04, and %l = 1-03. If Galactic ab- 
sorption is treated as a free parameter, an unrealistically low 
value is determined, and kT = lO^lo' 4 , 2 keV. In table[2]we 
list kT for each observation and the value from the com- 
bined dataset with Galactic absorption either fixed at the 
observed value, or treated as a free parameter. 

In the following subsections we present the X-ray de- 
rived profiles of the metal abundance, surface brightness, 
and temperature. 

3.1 Heavy element abundances 

Before assessing the observed temperature profile it is im- 
portant to determine the gas metal abundance. In principle, 
these two quantities may be decoupled with high quality 
data; however, the limited spectral resolution and signal-to- 
noise ratio means that the abundance is generally hard to 
constrain independently with radius, and one must adopt a 
mean value for the cluster. Vikhlinin et al. (2005) took ad- 
vantage of the superior spatial resolution of Chandra to show 




channel energy [keV] 



Figure 1. We show (upper panel) the observed X-ray spectrum 
of A1689 from the combined observations ID 540, 1663 and 5004 
in the 0.5 — 7 keV band, reduced from a 3' aperture. The fit (solid 
curve) is to a MEKAL model with Galactic absorption fixed at 
the observed value. We also show (lower panel) the residuals of 
the fit. 

that for a sample of nearby relaxed clusters there is a metal- 
licity gradient such that the abundance increases toward the 
center. However, the high roughly solar abundance level is 
actually in a region coinciding with the central galaxy. On 
the theoretical side, Arieli, Rephaeli & Norman (in prepa- 
ration) have performed a high resolution hydrodynamical 
simulation which uses a new approach to incorporate feed- 
back from galaxies on the intracluster gas, and have shown 
that including physical processes such as galactic winds and 
gas stripping yields a flat metallicity profile out to large 
radii (~ 600 h^ kpc). Based on XMM data, AM04 showed 
that there is no abundance gradient in A1689. Our deduced 
abundance gradient is shown in figure [2] There is a system- 
atic trend towards a higher value of A in the Chandra data 
(in accord with the result of XW02) than the corresponding 
XMM value, but this difference is not very significant from a 
statistical point of view, since most of the data points agree 
within the lrr range. The two sets of observations have the 
same exposure time, so the XMM data are more precise (due 
to the larger effective area). Our analysis indicates a fairly 
constant abundance over the region probed. In what follows 
we use the mean value of 0.4 solar for the central 3' region. 

3.2 Surface brightness and temperature analysis 

An HST/ACS image of A1689 overlaid on the X-ray map is 
shown in figure [3] the cD galaxy lies at the X-ray centroid. 

The symmetry of the surface brightness distribution al- 
lowed us to construct an azimuthally averaged profile in ra- 
dial bins. We determined an azimuthally averaged flux F in 
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Table 2. Mean gas temperature over a 3' region for each observation, with Galactic absorption either fixed at the observed value (Dickey 
& Lockman 1990) or left as a free parameter. 





N G = 2- 10 20 cm" 2 




Nq is a free parameter 




Obs. ID 


xl 


T [keV] 


xl 


T [keV] 


540 


1.16 


q 7+0.8 
y -'-0.5 


1.16 


lu -0.4 


1663 


1.1 


q+0.5 
3 -0.5 


1.01 


10 2+ ' 4 
±u - z -0.6 


5004 


1.19 


q n+0.5 
a ' J -0.4 


1.13 


10.50.gU 


combined 


1.03 


g 4 +o.4 


0.95 


10 7 +0 ' 3 
lU./_ 2 



0.9 
0.8 
0.7 



0.6 




Figure 2. The abundance profile as derived from the spectral 
deprojection fitting. We compare our Chandra analysis (circles), 
which gives a mean value of 0.42 (horizontal solid line), to the 
XMM analysis of AM04 (triangles) which gave a mean value of 
0.27 (horizontal dashed line). 



each radial bin by fitting the model to the data in each an- 
nulus and extracting the flux, from which the surface bright- 
ness S was found by dividing by the angular area of the bin. 
The error in S was calculated from that in the normaliza- 
tion parameter of the fitted model. Note that our profile 
of the surface brightness is ~ 4 times higher than that ob- 
tained by XW02 (see their Figure 6), but is consistent with 
that obtained by Mohr, Mathiesen & Evrard (1999, here- 
after MME99; see their Figure 13). 

The full extraction region consists of 50 annuli, each 
// 

of width 4 , plus nine annuli with varying widths (to keep 

// 

adequate S/N) beyond 200 , all centered on the cluster X- 
ray center (section 2.1). The maximum aperture of the data 
reduction region was limited by the outer edges of CCDs in 
practically all three observations. We removed the innermost 
ring because it corresponds to the inner region of the cD 
galaxy (see also Vikhlinin et al. (2005); e.g., the abundance 
in it was measured to be 0.9, contrary to our assumption of 
a constant abundance in each annulus. 

The reduction of the temperature profile was done in a 
similar way, but with a smaller number of annuli due to the 
lower signal-to-noise ratio. 




Figure 3. HST/ACS image of A1689 overlaid on smoothed color 
rendering of the X-ray emission (obs. ID 5004). The cD galaxy is 
seen to coincide with the X-ray centroid. The scale in the image 
is 0.5' ~ 65 h" 1 kpc. 

4 DERIVATION OF THE GAS AND MASS 
PROFILES 

4.1 Methodology 

In this section we present our procedure for studying the 
structure of the cluster using all the available data sets 
in combination. Since we have both lensing and surface 
brightness profiles, we do not need to assume any partic- 
ular parametrization or fitting formula, such as the com- 
monly adopted NFW mass profile or the j3 model for the 
gas density profile. We employ a model-independent ap- 
proach that is limited only by the resolution and accuracy 
of the data. As mentioned in the introduction, the cluster 
has only a small ellipticity, so the assumption of spherical 
symmetry is reasonable. This basic assumption allows us 
to reconstruct three-dimensional profiles from the observed 
two-dimensional ones. The lensing data yields in this way 
the density profile of the total mass (dark matter plus gas) , 
while the surface brightness profile depends on a combina- 
tion of the gas density and temperature. Thus, with one 
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additional relation, we can reconstruct the full gas and dark 
matter density profiles. We derive this additional relation by 
assuming hydrostatic equilibrium, consistent with the indi- 
cations discussed in section 1 that A1689 is well relaxed. 
Hydrostatic equilibrium involves the gravitational force and 
thus inserts a dependence on the total mass profile that cou- 
ples the constraints from lensing and from the X-ray data. 

We determined the best-fit values of our free parame- 
ters by fitting the lensing and X-ray surface brightness data 
simultaneously, ft is important to note that we did not in- 
clude the temperature profile data for reasons that are ex- 
plained later in this section. In our model-independent ap- 
proach, our free parameters are the values of the 3D profile 
of the total mass density and the gas mass density at sev- 
eral fixed radii (logarithmically spaced). The radial ranges 
of the free parameters of the total mass density and the gas 
density were determined by the range of the lensing data 
and the X-ray surface brightness data, respectively. Within 
these ranges the density values were linearly interpolated 
(in log-log) in-between the fixed radii. Beyond the last data 
point, each profile (i.e., of the total mass density and the 
gas density) was extrapolated. The total mass density pro- 
file was extrapolated as p(r) oc r -3 , in accordance with the 
asymptotic behavior of the NFW profile and also close to our 
best-fit core plus power-law model (section 5.2 below). The 
gas density was also extrapolated as p(r) oc r -3 , in order 
to have a constant / gas at large radii; this is also consistent 
with the previously deduced index of 2.9 (Xue & Wu 2000). 

In spherical symmetry, the hydrostatic equilibrium 
equation for an ideal gas can be integrated from a given 
three-dimensional radius r out to infinity (or some maxi- 
mum, cutoff radius). This yields the relation 



[p gas (r)T(r)] 



GMtotj^ r')prn pPgllB (r') , 
kr' 2 



(1) 



where p gas is the gas mass density, T is the gas temperature, 
M to t (^ r) is the total mass within r , p is the mean molecular 
weight, and m p is the proton mass. Of course, p depends on 
the abundances; the typical metal abundance is « | (which 
is expected when ejecta are from Type la supernovae winds 
(Silk 2003)). In A1689, A = 0.2- 0.49 was found by XW02. 

For given model parameters, the comparison to the lens- 
ing data was performed by first projecting the 3D profile of 
the total mass density using Abel integration 



k(R) = 



^crit 



Ptot{r)rdr 



R 2 



(2) 



where p to t (r) is the three dimensional total mass density and 
Sciit = t~c n "n i s the critical density for lensing, writ- 



4tvG D iD l3 

ten in terms of angular diameter distances D os (observer- 
source), D i (observer-lens), and D\ s (lens-source). Simul- 
taneously, we calculated the mass profile using 



M tot «r)=47r / p tot {r')r' 2 dr' 



(•3) 



We then used the mass profile and the gas density profile in 
the hydrostatic equilibrium equation ([T]), obtaining the tem- 
perature profile. The temperature and gas density profiles 
were then used (with the assumed abundances) to determine 
the emissivity: 



e(r) = n e (r)n H (r)A(T(r)) 



(4) 



where A(T) is the cooling function, which was obtained 
by considering all the relevant physical processes in the 
9-10 5 -3-10 8 K range. The cooling function was calculated by 
MEKAL (consistent with our fitting of the observed spectra 
in section [2~Tj) . We use the usual definitions, n = p gas /((j,m p ), 
n e = p gas /(p e m p ), and n H = p ga s/(jJ,Hm p ), which for 
A — 0.4 solar yield p = 0.55, p e = 1.05, and /j,h = 1-30. 

Having thus obtained the emissivity, we re-projected 
it using the Abel integral and obtained the X-ray surface 
brightness 



S(R) = 1.191 ■ 10" 



(1 + *) 



T 



e(r)rdr 



R 2 



(5) 



where e is the emissivity in erg s 1 cm , S is the X-ray sur- 



face brightness in erg s 
the units from erg s -1 cm 

D 2 

multiplying by 



cm 

2 



arcsec , and we converted 



to erg s 1 cm 2 arcsec 



2 by 
2 (where Dl 



4^dJ' usin S °l/Da = {l + z 
and Da are the angular diameter and luminosity distance, 
respectively), and upon conversion of radians -2 to arcsec -2 . 
The mean surface brightness was calculated in each annulus- 
shaped bin within the surface brightness data. The mean 
surface brightness Shin in a bin with inner and outer radii 
-Rin and -R ou t, respectively, is 

"-Rout 



Shin = 



{Rout Rin) 



RS(R)dR 



(6) 



Rin 



Finally, we compared this to the surface brightness data. 

Having calculated the observed quantities for each pos- 
sible set of model parameters, we then compared to the data 
using a \ 2 measure. The total x 2 is the sum of two terms, 
one from the lensing data and one from the X-ray surface 
brightness data: 



X 



2 . 2 
Xk+XS 



where 



K,(Ri) — Kdata(Ri) 



AKdata(Ri) 



and 



2 

Xs 



r Shin(Rj) — Sdata{Rj)~' 



E 



ASdata(Rj 



(7) 



(8) 



(9) 



Here x\ is the x 2 of the lensing data, obtained by comparing 
the model k of equation ([2]) with the N K data points Hdata 
and their errors Andata', similarly, xl is the x 2 of the X-ray 
surface brightness data, obtained by comparing the model 
S of equation (JSJ) with the Ns data points Sdata and their 
errors ASdata- The best-fit values of the parameters, which 
are the 3D profile of the total mass density and the 3D gas 
density profile, were obtained by minimizing the total \ 2 ■ 
From these two profiles we then derived the three dimen- 
sional temperature profile using equations (fTJ and ©. In 
order to perform an independent consistency check with the 
two-dimensional temperature data (which were not used in 
the analysis), we projected the 3D temperature, weighting 
it by the emissivity for comparison with the observed tem- 
peratures: 



r 2D (i?) 



lRnr)e{r) Tr 



-R? 



. dr 



lR<r) 



^Jr 2 -R? 



dr 



(10) 
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Weighting by the emission measure EM (which is propor- 
tional to Pg as ) gave similar results. 

The EM weighted temperature (Tem) and the emissiv- 
ity weighted temperature (Te) may not precisely reflect the 
actual spectral properties of the observed source. Vikhlinin 
(2005) states that the "spectroscopic" temperature is gen- 
erally lower than either of these temperatures, although 
his suggested correction formula for the effective spectro- 
scopic temperature of a multi-component thermal plasma 
yielded similar values to Tem and Te- Mathiesen & Evrard 
(2001) claimed that the observed spectroscopic temperature 
may be biased low due to the excess of soft X-ray emission 
from the small clumps of cool gas that continuously merge 
into the intracluster medium. Furthermore, Mazzotta et al. 
(2004) showed that from a purely analytical point of view 
the spectrum of a multi-temperature thermal model cannot 
be accurately reproduced by any single-temperature ther- 
mal model. It has also been found in simulations that the 
spectral temperature is in general smaller than Te by a fac- 
tor of 0.7 — 0.8 (Rasia et al. 2005). We note that total mass 
estimates obtained based solely on X-ray data together with 
the assumption of hydrostatic equilibrium depend strongly 
on the temperature. In our analysis, however, we do not use 
the temperature data for the fit but only for a consistency 
check. Therefore, our values for the total mass density are 
determined mainly by the lensing data. 

Our analysis focuses on the free-parameter method that 
does not assume a particular shape for the profile. In addi- 
tion, we used our analysis method to test the viability of 
either an NFW profile or a cored profile. 



4.2 The entropy profile 

It is interesting to evaluate the entropy and the adiabatic 
index profiles in order to probe the processes that govern the 
thermal state of intracluster gas. For a monoatomic ideal gas 
s = |fcln_ftT + so, where s — -S is the entropy per particle 
and K = Tn~ 2/:i . In this case A" is a combination of n 
and T that is invariant under adiabatic processes in the gas. 
More generally, T oc n r_1 under adiabatic changes, where 
r is the adiabatic index. In this paper we refer to K as the 
"entropy" . 

The gas at ~ 0.1i? v ir lies outside of the cooling re- 
gion, and is still sufficiently close to the cluster center where 
shock heating is minimized (Lloyd-Davis, Ponman, & Can- 
non 2000). Thus at this radius, it is less sensitive to specific 
models and assumptions. In the outer regions, r > 0.1-R v ir, 
the entropy is entirely due to gravitational processes, and 
the entropy profile is expected to be a featureless power law 
approaching K oc r (Tozzi & Norman 2001; Voit, Kay, 
& Bryan 2005). In the inner region the gas entropy pro- 
file is flattened (for ASCA and ROSAT data, see Lloyd- 
Davis, Ponman, & Cannon 2000; Ponman, Sanderson, & 
Finoguenov 2003; for various cosmological simulations see 
Voit, Kay, & Bryan 2005; for a theoretical analysis of joint 
X-ray and SZ observations, see Cavaliere, Lapi, & Rephaeli 
2005) . 



4.3 The effect of the cD galaxy 

A massive cD galaxy dominates the inner region of A 1689. 
Thus, it is important to check the effect of the cD galaxy on 
the cluster, and especially on the 2D temperature observed 
in the central region. Its mass profile can be analytically ex- 
pressed as follows. The mass profile is M c d = A/ Star s + Mdm, 
where the gas mass is neglected since it is dominated by the 
stellar mass at small radii and by the DM mass at large 
radii. We assume an NFW profile for the dark matter, and 

2 

express the stellar mass profile as M sta rs — Af sta rs o ■ , T v i , 

' L 1 r a ) 

where M s t ars ,o is the mass of the stars in the cD galaxy and 
a is a characteristic length-scale (Hernquist 1990). Typical 
values for the effective radius, R e , with a — R e /1.82, are 
R e ps 3 — 5 h -1 kpc. E.g., Jannuzi, Yanny, & Impey (1997) 
found an average R e = 4.6 h _1 kpc for five elliptic galax- 
ies; Brown & Bregman (2001) found an average R e = 3.2 
h _1 kpc for four elliptical galaxies (but see Arnalte Mur, 
Ellis, & Colless (2006) who found R e = 24.7 h" 1 kpc). Typ- 
ical values for the characteristic radius of the NFW profile 
of the DM (see also section 5.2 below) are 5 — 25 h" 1 kpc 
(Romanowsky et al. 2003). 

X-ray emission from the cD galaxy may be produced by 
the stars and stellar remnants, or by the hot gas. For the 
stellar component, the spectrum (dominated by X-ray bina- 
ries) is expected to be hard like that of the cluster (Tx ~ 10 
keV) and have a minor effect on the measured 2D tempera- 
ture, so we did not include it in our model. We considered 
the gas emission, which should give a soft spectrum (Tx ~ 1 
keV) (Brown & Bregman 2001). We adopted a double beta 
model for the gas density profile, as Matsushita et al. (2002) 
fitted to M87, where we redid the fit to the M87 data, al- 
lowing the parameters to vary freely except that we fixed 
the j3 of the extended component to be 0.47 (based on data 
from the ROSAT all sky survey: Bohringer et al. 1994). For 
a given gas mass density and total mass profile of the cD, 
we derived the temperature using equation Q. From the 
gas density and the temperature we obtained the emissivity 
using the MEKAL code, assuming a solar metal abundance. 



5 RESULTS 

In this section we present the results of our model- 
independent analysis of the lensing and X-ray data. The 
numbers of data points used in the analysis were N K = 26 
and N a = 58 for lensing and surface brightness, respectively. 
Accordingly, the ratio of the number of free parameters for 
the total mass to that for the gas mass profile was set to 
i/j ~ 26/58. As our standard case, we chose 6 free parame- 
ters for the total mass density and 14 free parameters for the 
gas mass density. These relatively small numbers of model 
points ensured fairly smooth profiles that effectively average 
over the noise in the data. The fit to the lensing and surface 
brightness data can be seen in fig.[l4]and fig. 1151 respectively. 
Both of the fits are very good, as expressed in the low re- 
duced x 2 , i.e., Xr = X 2 /dof. The reduced \ 2 of the fit to the 
data was 28.1/64. The contribution of each data set within 
the total, simultaneous fit to both, was Xr(lensing) = 4.7/20 
and Xr( s urface brightness) = 23.5/44. This shows that we 
achieved a good fit to both data sets, and thus did not need a 
larger number of model points. The low \t > especially in the 
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10' 

Radius [h~ 1 kpc] 

Figure 4. Results for the total mass density profile when different 
numbers of free parameters are used for the profiles of total mass 
(i) and gas (j). We consider i = 6, j = 14 (circles), i = 5, j = 11 
(triangles), and i = 4, j = 9 (squares). All three cases had the 
same, fixed values for the smallest and largest radius; these sets 
of points have been offset for display purposes, the triangles to 
the left and the squares to the right. 



10' 

Radius [h~ 1 kpc] 



Figure 5. Results for the gas mass density profile when different 
numbers of free parameters are used. Sets of points arc as in 
figure [4] The first and last points are offset as in figure [4] 



lensing case, suggests errors that are either overestimated or 
significantly correlated among the various data points, pro- 
ducing smoother data than expected given the errors. For a 
consistency check, we repeated the fits with different num- 
bers of parameters, namely i = 5,j = 11 and i = 4, j = 9. 
In figures [4] - [6] we show the profiles of the total mass den- 
sity, the gas mass density, and the temperature for the three 
free parameters sets, (i — 5,j = 11), (i = 6, j = 14) and 
(i = 4, j = 9). The figures clearly show that our results are 
insensitive to the precise number of adopted free parame- 
ters. To assess the degree of correlation between values of 
the projected fit parameters we have calculated the correla- 
tion matrix for i = 4,j = 9, which is specified in table lATI 
in the Appendix. 

We tested the effect of the form of extrapolation to 




10' 

Radius [h~ 1 kpc] 



Figure 6. Results for the 3D temperature profile when different 
numbers of free parameters are used. Sets of points are as in 
figure [J] The first and last set of points are again offset as in 
figure [4] 



large radii on the values obtained for the total mass den- 
sity and the gas mass density. We checked the results for 
large changes in the extrapolation power-law index, i.e. —2 
and —4 (compared to our standard assumption of p oc r~ 3 ). 
We used our standard number of free parameters, 6 and 14 
for the total mass density and for the gas mass density, re- 
spectively. In each case the same extrapolation index was 
used for the total mass and the gas mass density in order to 
ensure a constant gas fraction at large radii. We expected 
a large change in the free parameter at the largest radius, 
which is the closest to the extrapolation. Using an extrapo- 
lation index —2 gave a value which is lower by 21% and 11% 
for the total mass density and gas mass density, respectively, 
at this last point. For an extrapolation index —4 there was 
< 1% change in the total mass density value and 8% in the 
gas mass density. In the gas mass density there was also a 
change in the innermost radius, 29% and 36% for an extrap- 
olation index of —2 and —4, respectively; this change is not 
very significant due to the large errors at this radius. Taking 
a smaller number of free parameters lowers the errors of the 
gas mass density at the smaller radius and as a consequence 
also weakens the dependence on the extrapolation index. 
Thus, our results in general are not strongly dependent on 
the extrapolation index. 



5.1 Total mass density 

The values of our six free parameters of the total mass den- 
sity and the deduced 3D mass are shown in table ©. As 
mentioned above the X-ray temperature data were not used. 
The surface density data were used, but their effect on the 
total mass profile is limited since p ga s appears on both sides 
of equation ([1]). Thus, the lensing data basically determined 
the derived 3D total mass density and the 3D mass profiles. 

We also compared the total mass density profile ob- 
tained by the model-independent method to the that ob- 
tained by fitting particular models. We tested the NFW 
profile, 
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Table 3. Values of 3D mass density; errors are 1-<t confidence. 



r [h 1 kpc] p total [10 25 h 2 gr/cm 3 ] 



[10 1 



10.1 


qn+ 91 
9u -90 


0.91 ±0.68 


27.7 


39 ±11 


7.3 ±2.5 


75.9 


10.5 ± 1 


48.3 ±3.9 


208 


2.4 ±0.4 


252 ± 15 


568 


0.13 ±0.07 


679 ± 82 


1554 


0022+ 0049 
u.uuzz_ 0Q22 


956 ± 180 



Table 4. The values of the parameters of the two total mass 
density models, NFW and core. The errors are 1-cr confidence. 



Parameter 


NFW 


core 


po [1(T 25 h 2 gr/cm 3 ] 


9.6 ± 1.8 


24.4 ±4.5 


r s [h _1 kpc] 


175 ± 18 


91 ± 17 


n 




1.43 ±0.12 


X 2 /dof 


15.3/(26 - 2) 


13.3/(26 - 3) 



10 ,: 



» 


1 

i * t t 










/I 





10 10 10 

Radius [h~ 1 kpc] 

Figure 8. The derived 3D total mass profile. We compare our 
derived profile (dots) to the one derived by AM04 (squares) . They 
used an NFW profile based on data out to 693 h -1 kpc. Both 
profiles are shown with l-cx errors. The vertical dotted line is at 
0.1 r vir . 



(r/r s ) (1 + r/r s ) 2 
and a core model, 

P = (1 + (r°/r s ) 2 T ' (12) 

where r 3 is a scale radius. In the NFW profile po and r s were 
free parameters, and in the core profile po, r s , and n were 
free parameters. We chose these two profiles since they are 
frequently used and are significantly different both at small 
and at large radii. For simplicity, we fitted the two profiles 
only to the lensing data. 

A comparison of the profiles to the six values we ob- 
tained in our main analysis (figure [7] and table [4} yields a 
fairly close agreement and implies that neither of the profiles 




Radius [h 1 kpc] 



Figure 7. A comparison between the mass profile derived by the 
model-independent method to that using particular models. We 
compare the values obtained by our model-free analysis (points 
with error bars) to the results of assuming an NFW profile (solid 
curve) or a core plus power-law profile (dashed curve). 



is strongly excluded. Still, the measured profile as sampled 
by the six points is steeper at large radii than both models, 
in agreement with Broadhurst et al. (2005b). The core pro- 
file is too shallow at small radii and thus seems not to be a 
good fit to our deduced values of the mass density in this 
region. 

The NFW fit gave a concentration parameter Cjv = 
12.2]t5' 9 , where Cjv = r v i r /r s in terms of the virial radius 
r v ir and the characteristic radius r„ of the NFW profile. This 
is close to the value obtained by Broadhurst et al. (2005b), 
C N = 13.7±\-*, which was based only on a fit to the strong 
and weak lensing information. We also obtained from the 
NFW fit two characteristic radii for the cluster: r2oo, the 
radius inside which the average density is 200 times the crit- 
ical density, and r v i r , denned with a relative density of 116, 
which is the density expected theoretically at z = 0.183 in 
the ACDM model. The fit yielded r 2 oo = 1-71 h" 1 Mpc and 
r V ir = 2.14 h _1 Mpc. 

In figure [S] we plot our derived 3D total mass profile 
from eq. Q and the one derived in AM04 from XMM data. 
They used an NFW profile based on data out to 693 h _1 
kpc. The two profiles agree at the small radii but disagree 
at larger radii. Our highest-radius point agrees with their 
profile, but at this radius we had to extrapolate the AM04 
3D mass profile beyond their data range. 

5.2 2D mass profile 

The deduced 2D mass profile M2d(R) = 
2nT. cr it f Q R K(R')R'dR' , is compared with previous re- 
sults in Figure [9] These include results from strong lensing 
(Tyson & Fischer 1995) , the lensing magnification measured 
by the distortion of the background galaxy luminosity func- 
tion (Dye et al. 2001), the lensing magnification measured 
by the deficit of red background galaxies (Taylor et al. 
1998), the projected best-fit NFW model from X-ray data 
(AM04 via XMM), and the best-fit NFW model from weak 
gravitational shear analysis (Clowe & Schneider 2001; King 
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Figure 9. The derived 2D mass profile (6 diamonds) is compared 
to the gravitational lensing results from strong lensing (square, 
Tyson & Fischer 1995), distortion of background galaxy luminos- 
ity function (3 circles, Dye et al. 2001), deficit in number counts 
of red background galaxies (dashed curves, Taylor et al. 1998), 
projected best-fit NFW model from weak gravitational shear (tri- 
angles, King et al. 2002), and projected best-fit NFW model from 
X-ray data (solid curve, AM04) . 

et al. 2002). The latter analysis, although it is a weak 
lensing analysis, is the closest to the profile derived from 
X-ray data, but it appears to have underestimated the 
distortion signal derived in the other lensing analyses. This 
is probably caused by confusion between the cluster galaxies 
and the foreground or background galaxies. Thus, the X-ray 
analysis of AM04 (based on the X-ray temperature) gives a 
substantially lower 2D mass than the more reliable lensing 
analyses. This foreshadows the discrepancy that we find 
with the measured temperature (see section 5.5 below). 

5.3 Gas density profile 

The values of the derived gas mass density are shown in 
table [5] We compared our gas density profile to the one 
obtained by MME99. They fitted a single (3 model to the 
surface brightness profile; in clusters where they suspected 
the presence of cooling flows the fit was to a double /3 model 
(though with the same (3 in both elements). They used two 
types of clues for the presence of cooling flows. First, when 
fitting a single /3 model the surface brightness data points 
must be higher than the fitted curve in the inner region. 
Second, the cluster must appear relaxed, i.e., lacking obvi- 
ous asphericity or substructure that would indicate a recent 
merger. Judging by their successful fit to a single f3 model, 
they concluded that A1689 is not likely to have a cooling 
flow. Their fit is thus based on an assumed isothermal profile, 
and 0.3 solar abundance. Figure [10] shows excellent agree- 
ment between our obtained gas density profile and the one 
obtained by MME99. 

5.4 Gas mass fraction 

Having derived the total and gas mass densities, the gas 
fraction, / gas = M ga a/iV/totai, can then be determined upon 
integration. We can compare our value, f S s, B (r = 0.25r2oo) = 
0.0326 ± 0.0023 h" 3/2 , to Sanderson et al. (2003). They 
obtained for the 8 — 17 keV range (see their figure 6) 



Table 5 . Deduced values of the 3D gas mass density, and derived 
temperatures using the gas and total mass profiles together with 
the hydrostatic equilibrium equation. The errors are 1-cr confi- 
dence. 



r [h 1 kpc] 


Pgas [10~ 26 h - 5 gr/cm 3 ] 


T AD [keV] 


12.9 


o q+30.3 


10 1+ 30 ' 5 

— 10.1 


17.6 


7 9 +10.7 
'• z -7.2 


-i o 7+15.5 
lz - '-12.7 


24.0 


6.4 ±3.1 


13.0 ±4.7 


32.7 


5.06 ± 0.55 


14.0 ± 1.8 


44.5 


4.12 ±0.18 


14.1 ± 1.0 


60.7 


3.09 ±0.11 


14.7 ±0.75 


82.7 


1.96 ±0.055 


17.1 ±0.7 


113 


1.35 ±0.033 


17.8 ±0.7 


154 


0.97 ±0.023 


16.65 ±0.78 


209 


0.567 ±0.015 


16.94 ±0.96 


285 


0.35 ±0.012 


14.8 ± 1.0 


388 


0.19 ±0.0079 


13.1 ± 1.3 


529 


0.0982 ± 0.00067 


11.3 ±1.7 


721 


0.0586 ± 0.00057 


7.29 ±0.87 
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Figure 10. Our derived gas mass density profile (points with er- 
ror bars) compared to the profile derived by MME99 (solid curve). 

Table 6. Deduced values of the 3D gas mass at various radii. The 
errors are 1-cr confidence. 



t [h" 1 kpc] A/ gas [10 11 h" 5 / 2 M Q ] 



17.6 0.371^ 
24 0.71 ±0.57 

32.7 1.44 ±0.68 
44.5 2.93 ±0.64 
60.7 5.85 ±0.66 
82.7 10.93 ±0.65 
113 19.46 ±0.66 
154 34.63 ± 0.67 
209 59.3 ± 0.74 
285 96.83 ± 0.88 
388 151.5 ±1.3 
529 224.8 ± 2.5 
721 328.4 ±3.7 
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Figure 11. The gas mass profile derived from the best-fit gas 
density profile. All the values have error bars but some are too 
small to distinguish. These values are also listed in table itoT l, 



/ gas (V 0.25r 2 oo) = 0.041 b~ 3/2 . Our value of r 200 is 
slightly different since it was derived from the NFW fit 
to the lensing data, whereas their value was deduced from 
the X-ray data. They obtained r 2 oo = 2955 h^ = 2068.5 
h _1 kpc, and if we scale our result by their r 2 oo we get 
/ gas (r = 0.25r 2 oo) = 0.0359 ± 0.0024 h^ 3/2 . This value is 
still somewhat lower than their 8—17 keV value, but in 
good agreement if we look at the scatter of the gas frac- 
tion in their figure 4. In their figure 4 they plotted the val- 
ues of the gas fraction of different clusters at r = 0.3r 2 oo 
as a function of the cluster temperature. The gas frac- 
tion of hot clusters ranges from ~ 0.035 h" 3/2 to ~ 0.06 
h~ 3/2 , with some ~ 40% uncertainty. AM04 found for A1689 
/ gas (r = r 25 oo) = 0.044T0.0049 h" 3/2 , which is in agreement 
with our value of / gas (r = 7-2500) = 0.0385 ± 0.0026 h _3/2 . 
MME99 found for an ensemble of clusters with T > 5 keV 
a / g as(r = r 5 ao) = 0.075 h~ 3 ^ 2 , close to and slightly lower 
than the value in Sanderson et al. (2003), just as our gas 
fraction at r = 0.25r 2 oo was a little lower than the value 
derived by Sanderson et al. (2003) at that radius. 

Markevitch et al. (1999) examined two relaxed clusters, 
A2199 and A496, which are cooler than A1689 and have tem- 
peratures ~ 4 — 5 keV. They found / g as(»" = 0.5 h~ Mpc) = 
0.0569±0.0049 h~ 3/2 and / gas (r = 0.5 h' 1 Mpc) = 0.0559± 
0.006 h" 3/2 for A2199 and A496, respectively. These values 
are slightly higher than in other papers. They claimed that 
their values are consistent with others including MME99, 
who found that (/ gaa ) = 0.0566 h~ 3/2 for clusters with T < 5 
keV. However, note that the Markevitch et al. values are at 
riooo while the MME99 values are at rsoo- 

The gas mass fraction was also determined from SZ 
measurements with BIMA and OVRO (Grego et al. 2001). 
The latter authors measured f ga s(r = 65") and used a 
scaling relation in order to obtain / gas (7* = t"50o). They 
found for a sample of 18 clusters (assuming Q m = 0.3, 
Q A = 0.7) / gas = O.OSllo.o?? h _1 . These values are higher 
than obtained by X-ray measurements. For A1689, Grego 
et al. (2001) obtained / gaa (r = r 50 o) = 0.098±g;g|| h _1 , 
or more directly / gas (r = 65") = 0.068±^ 2 3 ft" 1 . Our 
results yield / gas (r = r 50 o) = 0.0482 h" 3/2 and / gas (r = 



65" = 140 h' 1 kpc) = 0.0222 ± 0.0009 h _3/2 , much lower 
than in their paper. Some systematic discrepancy in the 
gas fraction from X-ray and SZ measurements is expected 
(Hallman et al. 2007). The difference can be due to the 
fact that M gas (SZ) oc T" 1 (where T is the gas tempera- 
ture as deduced from spectral X-ray measurements; Grego 
et al. 2001), since (the thermal component of) the SZ ef- 
fect depends on the product of gas density and tempera- 
ture. Thus, an underestimation of the temperature (see sec- 
tion !5.5|l results in an M gas (SZ) that is higher than A/ gas (X), 
making / gas (SZ) higher than / gas (X). Expressing the dif- 
ference in terms of the gas mass, Grego et al. (2001) ob- 
tained Af gas (r = 65") = 4.61?;? ■ 10 12 h -2 M Q and we 
M gas (r = 65") = 2.9 • 10 12 h" 1 M Q . The ratio between 
the gas fraction derived by X-rays and the one derived by 
SZ also depends on dumpiness, where we define the clump- 
ing parameter as C = (n 2 ) 1//2 /(n e ). The gas mass measured 
via X-ray is M gas (X) oc C _1 , i.e., the actual gas mass is 
lower than that apparently observed if the cluster is clumpy 
but this is not corrected for. The gas mass measured via 
SZ is not proportional to the clumping. Thus, correcting 
for clumping would further decrease Af gas (X)/M gas (SZ) by 
a factor of ~ C. This may indicate that the level of clumping 
in this cluster is low. 

Finally, we also compared our results for the gas frac- 
tion to the simple theoretical models. In figure [12] we show 
the / gas obtained by the parametrized models, NFW (solid 
curve) or core (dashed curve) for the total mass density and 
(in both cases) a double /3 model for the gas mass den- 
sity, compared to the profile obtained from our model-free 
method (points with error bars). In general the models are 
not far from our reconstructed values. At very large radii, 
extrapolated beyond the data points, the gas fraction de- 
rived by the models would diverge from the real, suppos- 
edly constant value at large radii, since p gas oc r ~ 2,2 and 
Pdm /total r~ 3 . Instead, in our model-independent method 
the extrapolation of the gas density is p gas oc r -3 , in order 
to give an asymptotic / gas =const at large radii. At radii 
of ~ 200 — 500 h^ 1 kpc the / gas obtained by the model- 
independent method tends to be slightly lower than the one 
obtained by the parametrized models. At small radii the 
NFW model fits better than the core model, though the core 
gives a better overall fit to the lensing data (see table [4]|. 



5.5 2D temperature profile 

In figure [TH] we compare the 2D temperature obtained by 
our model-free method to the measured 2D temperature. 
The measured 2D temperature is substantially lower. It has 
been found in numerical simulations that the spectroscopic 
temperature is typically lower than the one determined by 
hydrostatic equilibrium by a factor of 0.7 — 0.8 (Rasia et 
al. 2005). Thus we also show in the figure the measured 2D 
temperature divided by 0.7 (upper solid curve). This fac- 
tor explains most of the discrepancy in A1689. Kawahara 
et al. (2006) used cosmological hydrodynamical simulations 
to explain this discrepancy. They found that local inhomo- 
geneities in the gas are largely responsible. 
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Figure 12. A comparison between the / gas profile derived by our 
model-independent method (points with error bars) to the pro- 
file derived using parametrized models. We consider two simple 
models, NFW (solid curve) and a core profile (dashed curve) for 
the total mass, in each case with a double beta profile for the gas 
mass. 



Figure 14. Our fit to the lensing data. The re profile obtained 
by our model-independent method (solid curve) is compared to 
the measured profile (points with error bars), where re is the 2D 
surface density in units of the critical density for lensing. 




150 200 250 
Radius [h~ 1 kpc] 

Figure 13. A comparison between the 2D temperature profile 
obtained from our model-independent method (dashed curve) to 
the profile measured from the X-ray spectrum (lower solid curve). 
Also shown are the measured temperature divided by 0.7 (upper 
solid curve) and the measured 2D temperature profile from AM04 
(circles with error bars). 



5.6 The entropy and polytropic index 

Previous theoretical work has shown that the entropy pro- 
files of non-radiative clusters approximately follow a power 
law with K(r) oc r 11 (Tozzi & Norman 2001; Bor gam et 
al. 2002; Voit et al. 2003). More recent detailed hydrody- 
namical simulations by Voit, Kay & Bryan (2005) confirmed 
the power law behavior of the entropy at r > 0.1r v i r . They 
also determined the normalization N, where K(r) / K200 = 
N(r/r 200 ) 1 - 1 , which was found to be 1.32 ±0.03 and 1.41 ± 
0.03 for SPH and AMR simulations, respectively. In the cen- 
tral region of the cluster, < 0.1r v ir, there are some claims 
for the existence of an entropy "floor", i.e., a limiting con- 
stant value. Voit et al. (2003) showed theoretically how an 
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Figure 15. Our fit to the surface brightness data. The X-ray 
surface brightness profile obtained by the model-independent 
method (solid curve) is compared to the measured profile (points 
with error bars). 



entropy floor can be attained; Voit, Kay & Bryan (2005) 
used simulations to find the level of this entropy floor. They 
found an entropy floor using an AMR code, but the value 
obtained with an SPH code was substantially lower (if at all 
present) . 

Observational comparisons of the entropy of clusters at 
small and at large radii show a departure from the predicted 
self-similarity. In the inner regions of low-mass clusters there 
appears to be a "floor" in entropy of ~ 135 keV cm 2 , with a 
higher value for high-mass clusters (Lloyd-Davies, Ponman, 
& Cannon 2000). This may represent a "preheated" mini- 
mum level which Kaiser (1991) speculated may be due to 
the effect of star formation in early galaxies which preheat 
the IGM through galactic winds. This idea is strengthened 
by the ubiquitous presence of gas outflows in observations 
of high-redshift galaxies (Franx et al. 1997, Frye & Broad- 
hurst 1998, Frye, Broadhurst & Benitez 2002, Pettini et al. 



A 1689 Mass and Gas Profiles 13 



20 




4 1 1 L , , , — . , . , — . . . 

10 1 10 z 10 3 

Radius [h~ 1 kpcj 

Figure 16. The 3D temperature profile as reconstructed by our 
model-independent method. 
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Figure 17. Our derived entropy profile (points with error bars) 
compared to the best- fit power law profile (solid curve). 

2001). Thus it seems natural to link this entropy floor with 
the winds from galaxy formation. 

In figure fTTl we plot K/K200, where K200 is the entropy 
at r — r 2 oo- We find r 20 o = 1-7 h" 1 Mpc and r 20 o = 0.8r v i r - 
The value of our derived entropy at 0.1r v i r agrees well with 
the one found by Lloyd-Davis, Ponman, & Cannon (2000), 
who used ROSAT and ASCA data. They assumed spherical 
symmetry and hydrostatic equilibrium, and a single beta 
model for the gas density along with a linear function for the 
temperature. For the DM density they used an NFW profile. 
Our derived entropy at 0.1r v j r , K(r — 0.1r v i r ) = 786 ± 33 
h -1//3 keV cm 2 , also agrees well with Ponman, Sanderson, & 
Finoguenov (2003), who also used ROSAT and ASCA data, 
assumed spherical symmetry and hydrostatic equilibrium, a 
single beta model for the gas density, and a linear function 
or a polytrope for the temperature. We fitted a power law 
to our entropy profile at all radii where we obtained the gas 
mass density (which includes points at r < O.lrvir). This 
fitted power law, which is shown in figure fTTl is K(r) / K200 = 
(0.95±0.05)(r/r 20 o)°' 82±0 ' 02 . 
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Figure 18. The adiabatic index deduced in this work (points with 
error bars), and the profile derived by Tozzi &; Norman (2001) 
(solid curve). 

Our X-ray data go up to ~ 0.34r v i r , so obtaining a fit to 
the (0.1 — l)r v ir range is not possible. In the available range 
~ (0.006 — 0.34)r v i r the power law is a close fit, though 
the index is somewhat lower than the theoretical value of 
1.1. Excluding the entropy points at small radii, r < 0.1r v i r , 
would not substantially change the value of the power index, 
but then there are only four data points at r > 0.1r v i r . In 
the inner part (r < 0.1r„i r ) if the entropy is flattened then 
this occurs only at r < Q.02r V ir- This is consistent with 
Tozzi & Norman (2001) who found that the flat entropy 
region m massive clusters (M ~ 10 15 M ) occurs only at 
r < 0.01 r 2 oo- Comparing our obtained values of the entropy 
at r < 0.1r v i r to those in Voit, Kay & Bryan (2005), they 
are closer to the values obtained in the SPH simulations 
than in the AMR simulations. Our derived normalization 
0.95 ± 0.05 is somewhat lower than obtained by Voit, Kay, 
& Bryan (2005) in both types of simulations. 

The deduced power law index is slightly lower than the 
value obtained from simulations (Tozzi & Norman 2001; 
Borgani et al. 2002; Voit et al. 2003; Voit, Kay, & Bryan 
2005). The disagreement can be explained by the fact that 
we did not assume an NFW profile. Assuming an NFW pro- 
file gives K(r)/K 20 o = (1.13 ± 0.39)(r/r 20 o) 1±0 2 . This is in 
good agreement with the "expected" power law index of 1.1. 
There is also a good agreement with the normalization ob- 
tained by simulations (Voit, Kay, & Bryan 2005). It has also 
been shown in simulations that as the amount of preheat- 
ing increases, the power law index becomes flatter at radii 
beyond the flat inner region (Borgani et al. 2005). 

In figure[l8]we plot the adiabatic index. The solid curve 
is taken from Tozzi & Norman (2001). They used an NFW 
profile for the DM and assumed no cooling. They also as- 
sumed an external, initial adiabat for the entropy. Our de- 
rived adiabatic index profile is in reasonable agreement with 
the profile derived by Tozzi & Norman (2001). 

5.7 The effect of the cD galaxy 

Figures [13] and [16] show that the 2D and 3D temperature 
profiles peak at r ~ O.lrsoo and decline at large radii. At 
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Table 7. The values of the best fitted parameters of the gas num- 
ber density of the cD assuming a double beta model for the gas 
mass density of cD. We fitted to the values taken from Matsushita 
et al. (2002) 



Parameter Value 



n e (l,0) [cm,- 3 ] 


117 +0.016 


r c> i [arcmin] 


S66+ 071 
u.doo_ 064 


0i 


439+0.833 


n e ,2(0) [cm -3 ] 


? ■ 10" 3 


r c ,2 [arcmin] 


5 41+0-31 
°- 41 -0.33 


#2 


0.47(fixed) 
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Figure 19. A comparison between the emissivity profiles of the 
cluster (dashed line) and of the cD galaxy (solid line). The cD 
emissivity profile is derived by assuming it is similar to that of 
M87. 



smaller radii the temperature profile also declines. This de- 
cline of the temperature profile at small radii occurs within 
the region that includes emission from the cD galaxy. Since 
the gas in the cD is denser and colder than the cluster gas, 
it might be the cause for this decline at small radii. In order 
to evaluate the effect of the cD on the cluster we modeled 
the emission of the cD. We estimated the gas mass density 
profile by taking the profile of M87, which is the cD in the 
nearest cluster of galaxies. We thus fitted a double j3 model 
(as Matsushita et al. 2002 did) to the surface brightness 
data in M87. The values of the best-fitted parameters of the 
gas number density of the cD are shown in table [7] Our ob- 
tained values agree well with the ones in Matsushita et al. 
(2002). For the evaluation of the cD mass profile we took 
R e = 7.8 h" 1 kpc, r s = 10 h" 1 kpc, M cD = 10 13 M Q , and 
Altars = 5 • 10 11 Mq (for the mass profile see section POl) . 

We checked the innermost part of the cluster r < 10 h _1 
kpc. This is out of the data range and requires an extrapo- 
lation, so we extrapolated using particular profiles. We used 
NFW for the DM profile and a double /3 model for the gas 
profile. In figure [19] we plot the emissivity of the cD galaxy 
and compare it to the emissivity of the cluster. The figure 
shows that the cD is dominant only at r < 4 h^ 1 kpc. Since 
the innermost surface brightness data point is at 4 h _1 kpc, 
the effect of the cD galaxy is likely minor, as expected in a 
rich cluster like A1689. 

We checked for an indication in the abundance profile 
that the cD contributes only at < 4 h _1 kpc (if at all). Note 
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Figure 20. A comparison of the abundance (in solar units) within 
4" annuli (points with error bars) to the average abundance of the 
cluster in a 3' aperture (solid line). The errors may be underesti- 
mated (see text). 

that the spatial resolution of Chandra is ~ 0.5", which for 
the redshift of A1689 is ~ 1.5 hp - * kpc. Also the cD galaxy 
coincides with the X-ray centroid within ~ 1.5" uncertain- 
ties (XW02). Figure |2"U1 indicates that the abundance may 
be higher, ~ 0.9, in the core than in the outer region. How- 
ever, this is not statistically significant due to the narrow 
annuli used which give high errors. The errors were calcu- 
lated here after freezing the temperature and the normal- 
ization. If these parameters are allowed to vary, this will of 
course increase the errors in the abundance. As expected 
from the low emissivity of the cD galaxy compared to that 
of the cluster, it is hard to see spectroscopically an indi- 
cation of the cD. We fitted a double temperature model, 
WABS (MEK AL+MEK AL) , to the 3' aperture, where the 
abundance of the first component (the cD galaxy) was fixed 
to solar and the second (the cluster) to 0.4 solar (the average 
value in the cluster). We also fixed the Galactic absorption 
to be 2- 10 -20 cm -2 as appropriate for this direction (Dickey 
& Lockman 1990). The resulting best fit is with two temper- 
atures, a cold one ~ 0.5 keV and a hot one ~ 10 keV. This, 
however, might be due to the real multi-temperature nature 
of the cluster gas. Extracting a smaller region, where the cD 
flux should be less diluted by the cluster flux, gave about the 
same two temperatures but with high errors. Thus, spectro- 
scopic fitting was inconclusive. 



6 DISCUSSION 

The increasing quality of X-ray and optical imaging data 
motivates renewed and more thorough examinations of the 
physical nature of galaxy clusters as revealed by the very dif- 
ferent processes of bremsstrahlung radiation and the gravi- 
tational deflection of light. In this paper we have examined 
the apparently relaxed cluster A1689, where only minimal 
substructure is evident from the dark matter, galaxy, and 
X-ray distributions and where the X-ray emission is smooth 
and symmetric. The longstanding claims of discrepancies in 
the total cluster mass estimated from these different kinds 
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of analysis can now be investigated with greater precision, 
fewer assumptions, and in a more model-independent way. 

There is a strong incentive to describe in analytical form 
the general nature of the total mass, gas, and temperature 
profiles. These quantities are inextricably bound up with the 
process of structure formation, including the nature of dark 
matter and the cooling history of the gas including interac- 
tion and merging of substructure. Thus, having an analyti- 
cal form of the profiles is very useful. The most commonly 
examined profile for the dark matter, based on N-body sim- 
ulations of collisionless dark matter, is the NFW profile. 
In contrast the most commonly fitted gas mass density or 
surface brightness profile, the j3 model, is essentially empir- 
ically based. The temperature profile is often derived from 
the analytical expression of the gas density using the poly- 
tropic equation of state, or the gas is taken to be isothermal 
since the temperature gradient was undetectable with the 
old X-ray satellites. This way of analyzing the cluster can 
lead to substantial errors in the DM density, gas density, 
and temperature values and the derived quantities such as 
the overall mass. 

As we discussed in the Introduction, the total mass pro- 
file can be independently determined from lensing measure- 
ments, or from X-ray measurements of the gas density and 
temperature profiles based on the assumption of hydrostatic 
equilibrium. (Clearly, the latter can be done also with spa- 
tially resolved S-Z measurements, when available.) The sec- 
ond basic assumption adopted in our analysis is spherical 
symmetry. Obviously, elongation along the line of sight is 
possible, but - as we have mentioned in the Introduction 
- this typically can introduce a 20% bias in the mass es- 
timate. For a detailed discussion of the impact of triaxial 
cluster morphology, see Gavazzi (2005). 

We have suggested here a model-independent method 
which uses free parameters and does not assume a specific 
profile for any mass component. The only hidden assumption 
is a linear (in log-log) interpolation between the free param- 
eters; this just assumes a reasonable degree of smoothness 
in the profiles. The profiles are also extrapolated beyond the 
data, but we showed that the results are insensitive to the 
detailed extrapolation. We specifically used a simple power 
law extrapolation for both the total mass density and for the 
gas mass density, with the same power law index assumed 
in order to approach a constant gas fraction at large radii. 
Our model-independent method is the best way to obtain 
the values of the important parameters of clusters at radii 
where there are good quality data. Our model-independent 
method can and will be applied to joint analyses of mea- 
surements of other clusters. 

Within the data range, the highest value of the 3D tem- 
perature (fig. 1161) is 2-3 times higher than the lowest value. 
Specifically, the temperature profile peaks at r ~ O.lrsoo and 
declines at both smaller and larger radii. The denser envi- 
ronment at the center of a cluster should naturally cause a 
decrease in the temperature there. We indeed find in this 
cluster that the 2D temperature decreases towards the cen- 
ter of the cluster, at r < 50 - 100 h' 1 kpc (figure . This 
can be due to several mechanisms. First, the high gas num- 
ber density in the center of the cluster causes a rapid loss of 
energy and a decrease in the temperature. As pressure sup- 
port decreases, the gas should gravitate towards the bottom 
of the cluster potential well, i.e., the center of the cluster, 



in a so-called "cooling flow" (Fabian 1994). A cooling flow 
should manifest in a detectable surface brightness enhance- 
ment in the X-ray emission from the central region of a clus- 
ter, ~ 100 — 200 kpc. Many studies have failed to find these 
cooling flows (MME99), at least not as expected in a simple 
or a standard cooling flow model (Fabian et al. 2001; Peter- 
son et al. 2001) . An explanation for this failure can be a more 
complicated mechanism which maybe includes reheating of 
the gas (e.g., by energetic particles; Rephaeli & Silk 1995), 
mixing, differential absorption, efficient conversion of cooling 
flow gas into low-mass stars, an inhomogeneous metallicity 
distribution (which is not the situation in A1689, at least 
not in terms of a radial gradient: fig. [2}, or disruption of 
cooling flows by a recent subcluster merger (see discussion 
in Fabian et al. 2001; Peterson et al. 2001). 

There are uncertainties about the criteria for finding 
cooling flows, and even with fixed criteria the decision 
whether cooling flows exist is still data dependent. At times 
the "solution" has been to exclude the inner, problematic 
region. Finding the value of the gas mass density in the 
model-independent approach bypasses these two problems, 
the unknown physical processes and the data dependence, 
as it allows for an analysis with fewer prior assumptions. 
Another way to reduce the core temperature is through the 
effect of the cD galaxy that is often anchored at the clus- 
ter center. The cD galaxy has a lower temperature and is 
denser than the cluster so it can in principle be the main 
cause for the low temperature at the center. However, based 
on our our initial, limited study it seems unlikely that the 
cD has a major effect in A1689, since the emissivity of the 
cD dominates that of the cluster only below r = 4 h _1 kpc 

(%■ OS). 

We have found a good fit to the observed 2D profiles 
of the lensing surface density and the X-ray surface bright- 
ness. Due to the smoothness of the profiles, we were able to 
fit the data with a relatively small number of free parame- 
ters. We have found a good agreement with previous results 
for all the parameters we have checked for A1689, including 
gas mass density, gas fraction profiles, measured tempera- 
ture, adiabatic index, and abundance (except compared to 
AM04). The total mass density profile we obtained was es- 
sentially determined directly by the lensing data alone, with 
< 1% differences introduced by using the X-ray data as well. 
This is the case since we did not use the X-ray temperature 
as part of the fit. 

We have shown that it is possible to obtain a model in- 
dependent 3D mass profile for which very good agreement is 
found between the lensing mass profile and the X-ray emis- 
sion profile. However, there is still a discrepancy between the 
temperature derived directly from (X-ray) measurements, 
and that deduced from a solution to the HEE, with the lat- 
ter ~ 30% higher at all radii, as has already been determined 
by Mazzotta et al.(2004) and Rasia et al. (2005). Other phe- 
nomena, such as bulk motions, turbulence, and nonthermal 
degree of freedom may contribute to the pressure, and thus 
reduce the temperature with respect to the one obtained 
from assuming that thermal gas pressure is the only con- 
tributor. Specifically, Faltenbacher et al. (2005) claim that 
based on their high-resolution cosmological simulations - in 
which they identified and analyzed eight clusters at z = - 
about 10% of the total pressure support may be contributed 
by random gas bulk motions, which may affect temperature 
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by up to 20%. It has also been suggested that X-ray lu- 
minous clumps of relatively low temperature may bias pro- 
jected temperature measurements downward (Kawahara et 
al. 2006). In any case, as we have already mentioned this 
temperature discrepancy is much smaller than the mass dis- 
crepancy by a factor of 2 — 4 encountered in previous (sep- 
arate X-ray and lensing) analyses. 

The derived entropy profile for A1689 has a power 
law form with no obvious central flattening, at least at 
r > 0.02r v i r , as expected for massive clusters from theo- 
retical work (Tozzi & Norman 2001). The existence and the 
level of the entropy floor is still not completely clear in simu- 
lations, since different simulation methods give different re- 
sults (Voit, Kay, & Bryan 2005). This is probably due to the 
limited resolution of simulations which is especially critical 
in the core of the cluster. Upcoming simulations with spa- 
tial resolution of ~ 5 kpc should give a better understanding 
of the entropy floor. Our derived radial slope is around 0.8 
rather than 1.1, flatter than the prediction. Since using sim- 
ple parametrized models for the DM and the gas gave a 
power law index of 1.1, we believe the value of 0.8 might be 
a more accurate result, which might be a reflection of a more 
complex gas dynamics and preheating history. Indeed, it has 
been suggested from simulations that preheating decreases 
the power law index in the region of our data (Borgani et 
al. 2005). 
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APPENDIX A: CORRLEATIONS BETWEEN 
PARAMETER VALUES 

In our model-independent approach we have determined pa- 
rameter values by fitting the projected mass, gas density, 
and temperature. Projection of the 3D quantities builds up 
correlations between their best-fit values in different radial 
bins. An illustration of the degree of such correlations is 
shown in table Al, which is the correlation matrix of the 
free parameters, taken from running four and nine free pa- 
rameters for the total and gas mass density respectively. The 
correlation in the deduced 2D temperature is shown in table 
A2. The elements of the 2D temperature correlation matrix 
correspond to the 12 values of the 2D radii (fig. 13). In both 
tables the correlations are quite strong between adjacent el- 
ements. In the 2D temperature correlation matrix there are 
also strong correlations between elements which are not ad- 
jacent but are at small radii, i.e. the upper left part of the 
table. 
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Table Al. The 13 X 13 correlation matrix for i = 4, j = 9 free parameters for the total mass and gas mass density, respectively. 





i=l 


i=2 


i=3 


i=4 


j=l 


j=2 


j=3 


j=4 


j=5 


j=6 


j=7 


j=8 


j=9 


i=l 


1.000 


-0.342 


0.197 


-0.072 


-0.026 


0.036 


0.034 


0.028 


0.040 


0.026 


0.026 


0.001 


-0.006 


i=2 


-0.342 


1.000 


-0.825 


0.273 


0.004 


-0.005 


-0.000 


0.010 


-0.008 


-0.051 


-0.074 


-0.005 


0.019 


i=3 


0.197 


-0.825 


1.000 


-0.599 


-0.002 


0.002 


-0.000 


-0.005 


0.030 


0.075 


0.091 


-0.001 


0.004 


i=4 


-0.072 


0.273 


-0.599 


1.000 


0.001 


-0.002 


-0.002 


-0.003 


-0.028 


-0.045 


-0.038 


0.018 


-0.074 


j=l 


-0.026 


0.004 


-0.002 


0.001 


1.000 


-0.749 


0.323 


-0.152 


0.057 


-0.026 


0.010 


-0.003 


0.002 


j=2 


0.036 


-0.005 


0.002 


-0.002 


-0.749 


1.000 


-0.613 


0.266 


-0.106 


0.047 


-0.018 


0.006 


-0.003 


j=3 


0.034 


-0.000 


-0.000 


-0.002 


0.323 


-0.613 


1.000 


-0.643 


0.238 


-0.107 


0.042 


-0.014 


0.007 


j=4 


0.028 


0.010 


-0.005 


-0.003 


-0.152 


0.266 


-0.643 


1.000 


-0.560 


0.234 


-0.095 


0.031 


-0.014 


j=5 


0.040 


-0.008 


0.030 


-0.028 


0.057 


-0.106 


0.238 


-0.560 


1.000 


-0.590 


0.223 


-0.075 


0.038 


j=6 


0.026 


-0.051 


0.075 


-0.045 


-0.026 


0.047 


-0.107 


0.234 


-0.590 


1.000 


-0.557 


0.174 


-0.086 


j=7 


0.026 


-0.074 


0.091 


-0.038 


0.010 


-0.018 


0.042 


-0.095 


0.223 


-0.557 


1.000 


-0.492 


0.226 


j=8 


0.001 


-0.005 


-0.001 


0.018 


-0.003 


0.006 


-0.014 


0.031 


-0.075 


0.174 


-0.492 


1.000 


-0.689 


j=9 


-0.006 


0.019 


0.004 


-0.074 


0.002 


-0.003 


0.007 


-0.014 


0.038 


-0.086 


0.226 


-0.689 


1.000 



Table A2. The 12 X 12 correlation matrix of the deduced 2D temperature for i = 4, j = 9 free parameters for the total mass and gas 
mass density, respectively. The raw and column numbers specify radial points (see fig. 13). 



1.0000 


0.9085 


-0.6717 


0.2874 


-0.1372 


-0.0423 


0.0604 


-0.0008 


-0.0233 


-0.0026 





0050 


-0.0005 


0.9085 


1.0000 


-0.3046 


0.1364 


-0.0854 


-0.0382 


0.0264 


0.0081 


-0.0008 





0044 





0100 


0.0093 


-0.6717 


-0.3046 


1.0000 


-0.3047 


0.1225 


0.0123 


-0.0773 


0.0206 


0.0515 





0162 





0111 


0.0219 


0.2874 


0.1364 


-0.3047 


1.0000 


-0.3307 


-0.1412 


0.1199 


0.0292 


-0.0119 





0156 





0363 


0.0304 


-0.1372 


-0.0854 


0.1225 


-0.3307 


1.0000 


0.2272 


-0.3406 


0.0566 


0.1892 





0501 





0213 


0.0626 


-0.0423 


-0.0382 


0.0123 


-0.1412 


0.2272 


1.0000 


0.8048 


0.0974 


-0.1890 





0403 





1515 


0.0908 


0.0604 


0.0264 


-0.0773 


0.1199 


-0.3406 


0.8048 


1.0000 


0.3025 


-0.0556 





0577 





1344 


0.1047 


-0.0008 


0.0081 


0.0206 


0.0292 


0.0566 


0.0974 


0.3025 


1.0000 


0.9125 





2147 





0351 


0.2234 


-0.0233 


-0.0008 


0.0515 


-0.0119 


0.1892 


-0.1890 


-0.0556 


0.9125 


1.0000 





4082 





1285 


0.2108 


-0.0026 


0.0044 


0.0162 


0.0156 


0.0501 


0.0403 


0.0577 


0.2147 


0.4082 


1 


0000 





6956 


0.2024 


0.0050 


0.0100 


0.0111 


0.0363 


0.0213 


0.1515 


0.1344 


0.0351 


0.1285 





6956 


1 


0000 


0.7934 


-0.0005 


0.0093 


0.0219 


0.0304 


0.0626 


0.0908 


0.1047 


0.2234 


0.2108 





2024 





7934 


1.0000 



